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The  paper  reviews  recent  studies  of  supersonic/hypersonic  flow  and  shock  wave  control  using 
plasma  energy  addition  and  magnetohydrodynamics  (MHD),  focusing  on  the  work  by  the  Princeton 
University  group  and  their  collaborators.  Applications  include  sonic  boom  mitigation  and 
optimization  of  scramjet  inlets.  Experimental  and  computational  studies  of  interaction  of  oblique 
shocks  with  laser-generated  blast  wave  and  thermal  wake  demonstrate  the  possibility  of  shock 
mitigation  in  a  repetitive-pulse  mode.  Analysis  of  interaction  parameter  for  MHD  control  of  cold 
hypersonic  flows  with  external  ionization  shows  that  significant  interaction  can  be  achieved  with 
energy-efficient  ionization  by  electron  beams.  The  maximum  achievable  interaction  parameter 
sharply  increases  with  increasing  Mach  number  and  altitude;  however,  interelectrode  arcing  may 
limit  the  performance.  For  MHD  control  of  scramjet  inlets,  nonequilibrium  electrical  conductivity 
is  created  by  electron  beams  injected  into  the  gas  along  magnetic  field  lines.  At  Mach  numbers 
higher  than  the  design  value,  the  shocks  that  would  otherwise  enter  the  inlet  can  be  moved  back  to 
the  cowl  lip  by  a  short  MHD  generator  at  the  first  compression  ramp.  To  increase  air  capture  at 
Mach  numbers  below  the  design  value,  a  heated  region  is  used  to  create  a  “virtual  cowl”  and  to 
deflect  flow  streamlines  into  the  inlet.  The  best  location  of  the  energy  addition  region  is  near  the 
intersection  of  the  nose  shock  of  the  vehicle  with  the  continuation  of  the  cowl  line,  and  slightly 
below  that  line.  Stretching  and  tilting  the  energy  addition  region  improves  performance.  By 
spending  only  a  few  percent  of  the  enthalpy  flux  into  the  inlet,  the  air  capture  and  engine  thrust  can 
be  increased  by  15-20%,  with  no  loss  in  specific  impulse. 
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Energy  addition  for  modification  of  the  external  flow  field  around  a  vehicle  is  a  well 
known  technique  and,  primarily,  has  been  studied  and  used  to  reduce  the  drag  on  objects  flying 
at  supersonic  speeds.  Historically,  energy  addition  in  the  supersonic  flow  was  proposed  as  an 

1  3 

alternative  technique  of  reducing  drag  on  a  flying  object.  ' 

In  principle,  with  a  proper  position,  shape,  and  power  of  the  thennal  source  relative  to  the 
external  surface  of  the  body,  one  can  substantially  decrease  its  aerodynamic  drag  and  thereby 
reduce  the  value  of  thrust  required  to  maintain  the  steady  flight.  Such  reasoning  was  discussed  in 
a  number  of  papers.1'8  The  problem  of  supersonic  flow  past  a  localized  region  of  heat  release 
was  treated  in  Refs.  4-7,  10-12. 

In  addition  to  drag  reduction,  other  applications  of  energy  addition  to  external  high-speed 
flows  have  been  suggested  recently.  Those  include,  for  example,  control  of  oblique  shock  waves 
and  optimization  of  scramjet  inlets,  control  of  shock  reflection  and  boundary  layer  separation, 
and  sonic  boom  mitigation. 

Most  recently,  Ref.  13  presented  an  overview  of  the  results  of  aerodynamic  flow  control 
by  energy  deposition  obtained  by  the  Princeton  group  and  their  collaborators  during  the  past  few 
years  and  insufficiently  reflected  in  the  recent  reviews  on  this  topic.11'  12  The  principal  topics 
included:13 

•  The  effect  of  a  plasma  (directed-energy)  air  spike  on  drag  reduction  and  effective  body 
geometry  in  hypersonic  flows. 

•  Modeling  of  “virtual  cowl,”  i.e.,  of  energy  addition  to  hypersonic  flow  off  the  vehicle  in 
order  to  increase  air  mass  capture  and  reduce  spillage  in  scramjet  inlets  at  Mach  numbers 
below  the  design  value. 

•  Numerical  studies  of  aerodynamic  forces  created  by  off-axis  heat  addition  upstream  of  a 
body.  These  forces  and  moments  can  be  used  for  supersonic/hypersonic  vehicle  steering. 

•  Experiments  and  numerical  modeling  of  pulsed  off-body  energy  addition  for  sonic  boom 
mitigation. 

In  the  present  paper,  we  briefly  review  just  two  of  the  listed  energy  addition  topics:  the 
“virtual  cowl”  concept,  and  sonic  boom  mitigation  by  energy  addition. 

The  issue  of  magnetohydrodynamic  (MHD)  control  of  hypersonic  flows  is  being 
extensively  studied  for  applications  in  hypersonic  airbreathing  propulsion.  Indeed,  the  geometry, 
size,  and  weight  of  scramjet-powered  hypersonic  vehicles  are  largely  dictated  by  the  need  to 
compress  the  ambient  low-density  air  upstream  of  the  combustor.14  The  optimum  geometry 
corresponds  to  the  well-known  shock-on-lip  (SOL)  condition:  the  compression  ramp  shocks 
converge  on  the  cowl  lip,  and  the  reflected  shock  impinges  on  the  inlet  shoulder.14  At  Mach 
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numbers  higher  than  the  design  value,  the  shocks  move  inside  the  inlet,  causing  local  hot  spots, 
loss  of  stagnation  pressure,  possible  boundary  layer  separation,  and  engine  unstart.  At  Mach 
numbers  lower  than  the  design  value,  a  portion  of  the  air  compressed  by  the  shock  misses  the 
inlet  (“spillage”),  resulting  in  thrust  decrease. 

To  avoid  perfonnance  penalties  at  off-design  Mach  numbers,  a  variable  geometry  inlet  can 
be  used.  An  alternative  approach  is  to  optimize  inlets  using  energy  addition  to  or  extraction  from 
the  flow.  Plasmas  and  magnetohydrodynamic  devices  may  offer  viable  optimization  schemes.15' 
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An  important  factor  in  MHD  control  is  that  in  the  flight  regime  of  interest,  static 
temperature  is  too  low  for  thermal  ionization  of  air.  Therefore,  nonequilibrium  ionization  in 
MHD  devices  must  be  applied.  ’  ’  The  energy  cost  of  ionization  imposes  very  rigid 

constraints  on  the  choice  of  ionization  methods:  only  the  most  energetically  efficient  ionization 
method,  that  by  electron  beams,  can  be  used  in  these  devices.  However,  even  with  the  most 
efficient  ionizer,  the  very  existence  of  ionization  cost  does  not  allow  the  ionization  degree  and 
the  conductivity  in  hypersonic  MHD  devices  to  be  very  high.  ’  '  ’  '  This  makes  the  task  of 
using  MHD  for  hypersonic  flow  control  and  power  generation  challenging  and  require  careful 
analysis  and  optimization. 

An  attractive  scenario  for  MHD  inlet  control  ’  ’  is  to  design  the  vehicle  for  a 

relatively  low  Mach  number  (for  example,  Mach  5  or  6),  and,  at  high  Mach  numbers,  move  the 
shocks  from  inside  the  inlet  back  to  the  cowl  lip  by  an  MHD  generator  device  placed  at  one  of 
the  compression  ramps.  The  advantages  of  this  method  are  the  following:  1)  the  vehicle  size  and 
weight  can  be  reduced  due  to  both  absence  of  a  mechanically  variable  inlet  and  the  shortening  of 
the  compression  ramp  (the  weight  savings  may  offset  the  weight  of  the  magnet  and  other  MHD- 
related  hardware);  2)  no  net  power  is  required  to  run  the  MHD  device,  since  the  generator  mode 
is  used,  and  the  power  requirements  for  air  ionization  can  be  minimized  with  electron  beams  as 
ionizers;  3)  flexibility  of  flow  control.  A  disadvantage  of  this  approach  is  that  the  inevitable 
Joule  heating  and  other  dissipative  processes  result  in  some  stagnation  pressure  losses. 

In  our  recent  work,  ’  '  ’  ‘  we  analyzed  inlet  control  at  Mach  numbers  higher  than  the 
design  value  by  means  of  an  on-ramp  MHD  generator.  The  principal  new  element  in  the  most 
recent  work  ’  is  that  MHD  regions  are  very  short  compared  with  both  length  and  width  of  the 
vehicle  forebody.  In  the  present  paper,  we  review  the  results  briefly. 

In  recent  papers,  we  also  suggested  and  analyzed  a  new  concept  of  increasing  the  mass 
capture,  dubbed  a  “virtual  cowl”.15' 35  The  essence  of  the  method  is  to  create  a  heated  region 
upstream  of  and  somewhat  below  the  cowl  lip.  The  incoming  flow  would  be  deflected  by  the 
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elevated-temperature  and/or  elevated  pressure  region,  causing  an  increased  mass  flow  into  the 
inlet.  The  heated  region  may  be  generated  by  supplying  microwave  or  RF  energy  to  a  volume 
pre-ionized  by  a  focused  laser  or  electron  beam.  Other  possible  means  include  hot  gas  or  plasma 
jets,  and  also  external  combustion.  Shooting  combustible  liquid  or  solid  pellets  (see  Ref.  21) 
upstream  may  be  an  interesting  option.  An  important  advantage  of  the  new  method  is  that  the  air 
entering  the  inlet  would  have  experienced  little  or  no  heating.  Thus,  irreversibilities  and 
stagnation  pressure  losses  associated  with  heating  can  be  minimized.  In  this  paper,  we  briefly 
review  the  results  of  computational  studies  of  the  virtual  cowl  concept,  and  present  recently 
perfonned  assessment  of  the  effects  on  engine  performance. 

2.  Energy  addition  in  supersonic  flows  for  sonic  boom  mitigation 

Shock  waves  generated  by  a  supersonic  object  are  known  to  produce  a  sonic  boom  on  the 
ground.  Sonic  boom  formation  has  been  the  main  obstacle  in  the  development  of  commercial 
supersonic  vehicles.  For  a  supersonic  aircraft,  the  near-field  shock  structure  is  a  complex  array  of 
shock  waves  that  originate  from  various  parts  of  the  aircraft.  In  the  far  field,  these  shock  waves 
coalesce  and  produce  an  N-shaped  pressure  signature  on  the  ground.  Although  the  sonic  boom  is 
a  far-field  phenomenon,  several  near-field  techniques  have  been  suggested  to  attenuate  the 
shocks  with  the  assumption  that  it  will  eliminate  or  partially  reduce  the  sonic  boom  on  the 
ground.  These  techniques  include  the  design  of  exceptionally  long  aircraft,  ’  use  of  underwing 
thermal  gradients,40  and  possibility  of  high  speed  oscillations  of  flight  velocity.41  Current  work 
focuses  mainly  on  the  development  of  long,  thin  aircraft  because  the  other  methods  seem 
unrealistic  due  to  practical  considerations.  Energy  addition  is  an  alternative  approach  to  sonic 
boom  mitigation.  It  has  been  argued  that  upstream  energy  addition  to  a  supersonic  vehicle  may 
weaken  the  shock  wave  and  may  prove  to  be  a  useful  approach  to  suppress  the  sonic  boom  on 
the  ground.  Recent  work  related  to  DARPA’s  Quiet  Supersonic  Platform  Program,  which  was 
conducted  at  Princeton,  addressed  these  issues  in  some  detail.  The  details  of  this  project  have 
been  presented  by  Miles  et.  al.  The  results  indicated  that  steady  state  off-body  energy  addition 
can  reduce  the  far-field  signature  primarily  by  suppressing  the  far-field  coalescence  of  the  shock 
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waves  originating  from  the  various  parts  of  the  vehicle. 

The  dynamic  effects  of  energy  addition  were  investigated  with  a  time-marching  Euler  code 
based  on  2nd  order  MacCormack  method  on  rectangular  grid.  Preliminary  modeling  was 
perfonned  in  two  dimensions  with  the  energy  added  along  a  line  across  the  two-dimensional 
domain  and  the  effects  seen  on  the  shock  structure  from  a  wedge.  The  simplicity  of  this 
configuration  allows  the  code  to  follow  the  dynamics  of  the  interaction  of  the  energy  addition 
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with  the  bow  shock.  An  example  of  the  interaction  of  a  series  of  energy  addition  pulses  is  shown 
in  Fig.  1.  These  are  five  frames  from  a  ‘movie’  that  show  the  dynamics  of  the  interaction.  Fig.  2 
shows  the  impact  of  those  pulses  on  the  pressure  measured  at  the  bottom  of  the  frames  shown  in 
Fig.  2  (z  =  -1.6m). 


Fig.  1.  Two-dimensional,  time-accurate  model  of  pulsed  energy  addition.  Flow  is  from  left  to 
right.  The  energy  addition  regime  modulates  and  weakens  the  bow  shock 


Fig.  2.  Maximum  relative  pressures  at  the  bottom  cross  section:  z=-1.6  m  on  Fig.  1.  Rectangles 
show  the  energy  addition  pulses  in  relative  units.  The  bow  shock  pressure  is  modulated  in  time 
(solid  line),  but  averaged  pressure  amplitude  from  the  bow  shock  is  20%  lower  than  steady-state 
one  without  the  heat  addition  (shown  by  the  arrow).  The  pressure  associated  with  the  shock 
generated  by  the  heat  addition  is  shown  by  the  dotted  line  and  maximal  pressure  from  the  heat 
source  is  significantly  lower  than  one  from  the  model.  Pulse  duration  r  ,  =  2.5  ms;  frequency 

f=125  Hz. 


In  the  absence  of  energy  addition,  the  bow  shock  produces  a  Pmax/P0  of  0.027  which  is  shown  by 
the  arrow  on  the  left  side  of  Fig.  2.  In  the  presence  of  the  pulsed  energy  addition,  the  Pmax/Po 
from  the  bow  shock  oscillates,  but  is  consistently  lower  in  value  than  it  is  in  the  absence  of  the 
energy  addition.  The  average  reduction  in  bow  shock  strength  is  about  20%.  The  energy  addition 
process  itself  generates  an  oscillating  shock,  but  that  is  significantly  weaker,  as  shown  by  the 
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dotted  line.  It  is  expected  that  the  addition  of  energy  by  pulses  such  as  this  will  reduce  the  total 
power  requirements  for  modifying  the  far-field  signature,  and  it  may  introduce  new  effects  that 
further  suppress  the  far-field  coalescence  and  reduce  sonic  boom.  These  dynamic  far-field 
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effects  cannot  yet  be  predicted  and  have  not  yet  been  measured. 

The  above  example  explicitly  demonstrates  the  significance  of  energy  addition  to  control 
flow  fields  around  objects  flying  at  supersonic  speeds.  Among  various  methods  of  energy 
addition,  laser  initiated  optical  discharges  have  been  found  to  be  an  efficient  way  to  control  the 
position  and  structures  of  shock  waves.11'13'  43  45  In  this  regard,  laser  based  energy  addition 
techniques  can  find  an  important  application  to  sonic  boom  mitigation.  To  demonstrate  this 
approach  and  to  investigate  the  dynamic  effects  of  energy  addition  upstream  of  the  shock  wave 
produced  by  a  model  in  the  flow,  experiments  were  performed  under  DARPA’s  Quiet 
Supersonic  Platfonn  program  at  Princeton. 

The  experimental  program  was  established  to  validate  model  predictions  for  the  dynamic 
energy  addition  and  to  examine  the  real  time  near-field  interaction  of  the  energy  addition  with 
the  shock  wave  produced  in  the  flow.  A  Mach  2.4  nozzle  was  designed  to  perform  experimental 
studies.  A  model  was  placed  in  the  test  section  [length=6",  width=2.15",  height=1.9"]  to  produce 
a  shock  wave.  The  static  pressure  in  the  test  section  was  near  atmospheric  to  pennit  localized 
energy  addition  by  laser  induced  breakdown.  A  10  Hz  Nd:YAG  laser  was  employed  to  add 
energy  into  the  flow.  The  pulse  energy  was  about  350  mJ/pulse,  and  a  100  mm  focal  length  was 
used  to  obtain  a  breakdown  in  the  region  of  interest  upstream  of  the  model.  A  Princeton 
Scientific  Instrument  ultra-fast  framing  camera  (PSI-4)  was  used  to  capture  the  shock  wave 
interactions  in  the  test  section  as  the  energy  is  added  to  the  flow.  Fig.  3  shows  the  experimental 
shadowgraph  images  captured  with  the  desired  flow  in  the  tunnel.  Flow  is  from  left  to  right,  and 
the  laser  was  focused  at  the  position  mentioned  in  Fig.  3.  The  frame  integration  time  was  2  ps. 
Frames  in  Fig.  3  clearly  show  the  thennal  spot  along  with  its  associated  shock  wave  propagating 
with  the  flow.  As  the  thermal  spot  passes  through  the  oblique  shock,  a  significant  distortion  and 
weakening  of  the  shock  can  be  seen  in  frames  4  to  7.  Full  details  on  these  experiments  can  be 
found  in  Ref.  43.  It  is  important  to  note  that  for  an  effective  sonic  boom  reduction,  the  position 
of  the  laser  breakdown  spot  should  be  optimized  so  that  the  shock  wave  generated  by  the  thennal 
spot  weakens  before  it  reaches  the  bow  shock.  In  this  case,  the  interaction  of  the  thermal  spot 
would  have  a  more  pronounced  impact  on  the  bow  shock. 
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(13)  (14) 

Fig.  3.  Shadowgraph  images  of  the  interaction  between  the  laser  thermal  spot  and  blast  wave  and 
the  oblique  shock  wave  from  a  wedge  in  a  Mach  2.4  flow.  Images  are  4  ps  apart. 

Modeling  was  conducted  at  the  flow  conditions  used  for  the  experimental  studies  (M=2.4, 
Po=l  atm,  To=136  K,  y=1.4).  A  time  marching  two  dimensional  Euler  code  was  used  for  this 
purpose.  The  energy  deposition  region  was  assumed  to  have  a  cylindrical  shape  and  a  Gaussian 
distribution  in  the  (x,y)  plane  with  an  effective  radius  of  0.125  cm.  This  assumption  was  based 
on  the  fact  that  a  laser  spark  produces  a  cylindrical  shock  wave  and  a  rarefied  thennal  spot  that 
both  move  downstream  with  the  flow. 

Fig.  4  presents  the  computational  results  describing  the  dynamics  of  interaction  of  the 
energy  addition  with  the  shock  wave  generated  by  the  model.  The  energy  addition  itself 
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generates  a  shock  wave  that  grows  as  it  propagates  in  the  flow.  Frames  in  Fig.  4  clearly  show  the 
complex  interaction  of  the  entropy  spot  with  the  oblique  shock  and  also  show  attenuation  of  the 
shock  wave  as  it  passes  through  the  entropy  spot.  Fig.  4  also  shows  the  way  the  laser  generated 
shock  wave  gets  reflected  from  the  model  surface,  which  makes  the  dynamics  of  interaction  even 
more  complex. 

Calculations  were  performed  for  the  interaction  of  the  flow-convected  thermal  spot  with 
and  without  the  accompanying  shock  wave  with  the  oblique  shock  from  the  model.  A  remarkable 
difference  between  the  two  cases  can  be  seen  in  terms  of  pressure  variation  along  the  surface  of 
the  test  section,  as  shown  in  Figure  5.  In  the  case  when  the  laser  spot  is  focused  10  mm 
upstream  and  10  mm  below  the  model,  the  pressure  variation  is  maximum.  The  rise  in  pressure  is 
due  to  the  interaction  of  the  shock  wave  surrounding  the  thennal  spot  with  the  oblique  shock 
wave,  whereas  the  pressure  decreases  when  the  thermal  spot  interacts  with  the  model  shock. 
Since  pressure  calculations  were  made  in  the  near  field,  only  2.5  cm  below  the  model,  the 
interaction  of  the  shock  waves  with  the  expansion  fan  located  at  the  far  right  edge  of  the  model 
affects  the  pressure  in  a  complex  fashion.  Due  to  this,  several  fluctuations  in  the  pressure  curve 
can  be  seen  in  Fig.  5.  As  the  laser  spark  is  moved  farther  away  from  the  model  edge  (40  mm 
upstream),  the  shock  wave  around  the  thennal  spot  becomes  weaker  and  the  conesponding 
pressure  variations  become  smaller  in  amplitude,  as  can  be  seen  in  Fig. 

5.  For  the  isobaric  case  where  the  pressure  inside  the  energy  deposition  region  has  been 
assumed  constant  (thennal  spot  without  the  surrounding  shock  wave),  only  the  interaction  of  the 
thermal  spot  has  been  analyzed,  and  the  pressure  associated  with  the  oblique  shock  wave,  unlike 
the  previous  two  cases,  was  reduced  without  any  large  fluctuations  (Fig.  5).  Thus,  the  position  of 
the  laser  breakdown  spot  should  play  a  vital  role  in  the  mitigation  of  shock  waves  and  of  the 
sonic  boom  from  a  supersonic  vehicle.  Far-field  experiments  will  be  required  to  demonstrate  the 
effectiveness  of  this  technique. 
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Fig.  4.  Two  dimensional,  time  accurate  model  of  laser  spark  and  the  oblique  shock  wave 
interaction  (density  contours).  Flow  is  from  left  to  right. 
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Fig.  5.  Top  -  computational  predictions  of  maximum  relative  pressures  at  the  bottom  of  the  test 
section  for  various  energy  addition  positions;  bottom  -  geometry  indicating  the  position  for 
localized  energy  addition  to  predict  pressure  variations  at  the  bottom  of  the  test  section 

The  results  of  both  experiment  and  numerical  modeling  clearly  demonstrate  strong 
influence  of  the  pulsed  power  deposition  in  supersonic  flow  on  shock  wave  intensity  in  the  near 
field  of  the  body.  But  whether  the  problem  of  sonic  boom  mitigation  in  the  far  field  can  be 
solved  by  such  means  requires  additional  theoretical  and  experimental  investigations.  It  is 
possible  that  the  ‘quiet’  power  deposition  (a  long,  thin  heated  region  with  low  power  input  per 
unit  volume)  suggested  by  Kolesnichenko46,47  and  by  Georgievskii  and  Levin48  can  be  useful  not 
only  for  drag  reduction,  but  also  for  sonic  boom  mitigation,  because  in  that  case  most  of  the 
power  is  spent  on  creation  of  the  thennal  wake,  rather  than  on  generating  shocks. 


1.  Basic  analysis  of  MHD  flow  control  in  cold  hypersonic  flows:  the  interaction 

29 

parameter  — 

To  estimate  MHD  effects  in  relatively  cold  hypersonic  flow  and  boundary  layer,  the 
momentum  equation  for  the  gas 

p  —  =  -Wp  +  /  x  B+pAu  (1) 

dt 

should  be  analyzed  jointly  with  other  appropriate  gas  dynamic  and  plasma  kinetic  equations,  and 
with  the  generalized  Ohm’s  law 
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y  =  (7|«+^x5  +  M(yxfi)xfi  (2) 

where  E"  =E+uxB  is  the  electric  field  in  the  reference  frame  moving  with  the  gas,  and  Q  and 
Q,  are  the  electron  and  ion  Hall  parameters: 


Q 


eB 

mkenn  ’ 


Q.  = 


eB 


M„k„ji 


(3) 


The  second  term  in  the  right-hand  side  of  Eq.  (2)  represents  the  Hall  effect,  and  the  third  term, 
nonlinear  with  respect  to  B ,  represents  the  ion  slip.  In  the  analysis,  we  will  disregard  negative 
ions  and  also  assume  quasineutrality,  ne  =  n+  if!  n  . 

The  MHD  interaction  parameter,  also  referred  to  as  the  Stuart  number,  S,  is  the  ratio  of 
ponderomotive  (Ampere)  and  inertia  forces.  With  ion  slip  correction, 


5  = 


aB2L 

(l  +  )pu 


(4) 


where  L  is  the  length  of  MHD  region.  When  ion  slip  is  small,  f )  Q  □  1 ,  the  Stuart  number 
increases  with  magnetic  field  as  B2.  However,  at  very  strong  magnetic  fields  and  in  low-density 
gases,  when  QQ  □  1 ,  the  interaction  parameter  reaches  its  asymptotic  value  independent  of  B: 
S->(L/u)kn+n+  (5) 

The  physical  meaning  of  Eq.  (5)  is  that  the  maximum  ion  momentum  change  occurs  when  they 
are  essentially  stopped  by  the  strong  transverse  magnetic  field,  and  then  momentum  transfer 
from  the  ions  to  the  gas  is  limited  by  the  number  of  collisions  of  a  neutral  molecule  with  ions. 
Indeed,  the  right-hand  side  of  Eq.  (5)  is  simply  the  ratio  of  the  flow  residence  time  in  the  MHD 
region  to  the  mean  time  for  a  molecule  to  collide  with  an  ion. 

As  discussed  in  our  earlier  work,  at  flight  Mach  numbers  up  to  about  Mach  12,  thermal 
ionization  cannot  provide  an  adequate  electrical  conductivity,  and  the  need  to  spend  power  on 
artificial  ionization  severely  limits  performance  of  hypersonic  MHD  devices.  ’  ’  ’  Indeed, 

the  work  done  on  an  electron  by  the  induced  Faraday  electric  field  during  the  electron’s  lifetime 
with  respect  to  dissociative  recombination  with  ions  is: 

1  (sEY  ( ekuBY 

s  =  eEVdr - =  y  ’  (6) 

P’h  mken(Jnne  mkmf3nne 


eE 

where  k  is  the  load  factor,  E  =  kuB  is  the  induced  electric  field,  Vedr  = -  is  the  electron  drift 

mkenn 

velocity,  and  p  is  the  dissociative  recombination  rate  constant.  For  efficient  self-powered 
operation  of  the  MHD  device,  s  must  be  substantially  larger  than  the  energy  cost,  W{ ,  of  a  newly 
produced  electron.  Specifically,  the  ratio 
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must  be  limited  to  a  number  less  than  1.  This  obviously  limits  the  electron  density  ne  and  the 
conductivity,  so  that  the  maximum  interaction  parameter  per  unit  length  is: 


S/L  = 
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e  kMn 

V  J 
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*4  (8) 


Assuming  that  the  MHD  device  is  located  downstream  of  a  single  oblique  shock  with  a 
small  turning  angle  6 ,  we  obtain  the  final  expression  for  the  maximum  Stuart  number  per  unit 
length: 
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where  q  is  the  flight  dynamic  pressure,  and 

f(M,0)  =  y— 

y  + 1 

32 

+ - r - 1  2  ,  =7 

(y  +  l)  (y  +  l)M Xand  +  +  l)  M‘tan‘^  +  16 

g(M,9)  =  1 - 

[y  +  P)M  \nn9  +  y^y  +  \)~  M1  tan2  9  + 16 

MHD  perfonnance  as  expressed  by  Eq.  (9)  is  inversely  proportional  to  the  energy  cost  of 
producing  an  electron,  .  Thus,  minimization  of  W,  is  critical.  As  shown  in  our  earlier  work,  ’ 
15,  17,  23  e|eclron  beams  represent  the  most  efficient  nonequilibrium  method  of  ionization,  with 
W,=  34  eV  ,  and  in  what  follows,  we  will  assume  that  ionization  is  done  by  electron  beams,  with 
W,  =  34  eV  ,  and  that  a  is  limited  to  0.3. 

The  maximum  interaction  parameter  (9)  increases  very  rapidly  with  increasing  magnetic 
field  and  flight  speed,  and  with  decreasing  dynamic  pressure.  The  very  sharp  m2  dependence  is 
due  to  both  increase  in  the  Faraday  e.m.f.  with  flow  velocity  and  to  the  decrease  in  gas  density  at 
constant  q.  Figure  6  shows  SI L  calculated  with  Eq.  (9)  versus  flight  Mach  number  at  4  different 
magnetic  field  strengths  and  2  values  of  flight  dynamic  pressure.  The  load  factor  value  was 
constant  at  AH). 5. 

Note  that  in  Faraday  MHD  devices,  Hall  current  is  eliminated  by  segmenting 
electrodes.49  However,  as  the  longitudinal  Hall  electric  field  increases,  arcing  between  the 


_ 32 _ 
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electrode  segments  can  occur.49  The  arcing  would  essentially  result  in  a  continuous-electrode 
Faraday  device,  with  dramatic  reduction  in  performance.  Denoting  the  value  of  Hall  field 
corresponding  to  the  arcing  threshold  as  Ec ,  the  upper  limit  for  the  interaction  parameter  per  unit 
length  becomes  with  this  constraint, 


Slim  ...  a Jli  SeEcU oo 

L  2  pWtq  [  2(1  ~k)kn+q^ 


(10) 


This  formula  does  not  explicitly  depend  on  B  field  and  has  a  weak  dependence  on  flight 
conditions,  because  at  each  flight  regime,  B  field  is  adjusted  to  satisfy  the  Ec  constraint. 

Fig.  6  shows  the  Hall  field-limited  interaction  parameter  of  Eq.  (10)  at  two  values  of  Ec. 
While  the  performance  is  quite  good  at  Ec  =  50  kV/m ,  it  is  unacceptably  low  with  Er  =  5  kV/m  . 
Therefore,  determining  the  threshold  field  for  intersegment  arcing  is  critical  for  assessing 
performance  of  hypersonic  MHD  devices. 


q=1  Atm;  9=10°;a=0.3 

1  -  B=1  T 

2  - B=2  T 

3  -  B=4  T 

4  -  B=7  T 
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Fig.  6.  Maximum  MHD  interaction  parameter  per  unit  length  versus  flight  Mach  number.  Solid 
lines  -  interaction  parameter  at  different  values  of  magnetic  field,  with  no  constraint  on  Hall 
field;  dashed  lines  -  upper  limit  of  interaction  parameter  imposed  by  constraint  on  maximum 
allowed  Hall  field. 


For  MHD  control  of  the  boundary  layer,  one  can  define  an  MHD  interaction  parameter 
with  respect  to  shear  stress  at  the  wall: 
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where  u'  is  the  friction  velocity  defined  by  pu'2  =  t  ,  where  the  wall  shear  stress  can  be  related  to 
freestream  conditions  p,  u  ,  and  the  friction  drag  coefficient 


r  = 


1  2 
—  pu  c, 
2  f 


(12) 


If,  for  estimates,  we  assume  that  the  conductivity  a  and  the  density  p  near  the  wall  are  equal  to 
those  in  the  core  flow  (this  can  be  ensured  by  wall  cooling  and  by  contouring  the  profiles  of 
ionizing  beams),  then: 


S  =  —  x 


<T  B2L 
pu 


=  —xS 


(13) 


In  the  range  Re*=105  -107,  the  MHD  interaction  parameter  with  respect  to  wall  shear  stress 
greatly  exceeds  S:  ST  ~  305 .  This  may  have  interesting  implications  for  wall  friction  and 
transition  control.  Of  course,  MHD  effects  on  the  boundary  layer  can  be  quite  complex.  If  the 
load  factor  k  is  not  close  to  1  or  0,  then  the  heating  effects  should  be  comparable  to  those  of 
ponderomotive  JxB  forces  both  in  the  core  flow  and  near  the  wall.  Also,  if  the  MHD  region 
extends  into  the  core  flow,  the  changes  in  the  core  flow  would  also  affect  the  boundary  layer. 


4.  MHD  inlet  control  at  M>MriPSiffn 

An  attractive  scenario  for  MHD  inlet  control  that  has  been  analyzed  by  several  groups,15' 
’  "  ’  ’  is  to  design  the  vehicle  for  a  relatively  low  Mach  number  (for  example,  Mach  5  or  6), 
and,  at  high  Mach  numbers,  move  the  shocks  from  inside  the  inlet  back  to  the  cowl  lip  by  an 
MHD  generator  device  placed  at  one  of  the  compression  ramps.  This  method,  while  providing 
some  flexibility  of  flow  control,  would  not  require  any  net  power  to  run  the  MHD  device,  since 
the  generator  mode  would  be  used,  and  the  power  requirements  for  air  ionization  can  be 
minimized  with  electron  beams  as  ionizers.  ’  '  ’  ’  Principal  disadvantages  of  the  method 
include  stagnation  pressure  losses  due  to  the  inevitable  Joule  heating,  and  also  weight  of  the 
magnet  and  other  hardware. 

In  our  recent  paper,  we  demonstrated  that  this  approach  to  MHD  inlet  control  can  yield 
an  acceptable  inlet  performance,  while  operating  in  the  mode  of  net  power  extraction.  We 
showed  that  the  MHD  region  should  be  quite  short  (typically  25-30  cm  along  the  flow)  and  be 
located  as  far  upstream  as  possible.  While  encouraging,  those  results  relied  upon  the 
assumption  that  the  magnetic  coil  is  very  large:  the  coil  diameter  was  assumed  to  be  equal  to  the 
forebody  width.  Minimizing  the  magnet  size  by  optimizing  the  MHD  region  would  be  important 
for  the  practicality  of  the  concept. 
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In  this  paper,  as  in  Ref.  37,  we  build  upon  our  earlier  work  and  attempt  to  optimize  the 
MHD  region,  so  that  the  magnet  size  is  minimized,  while  restoring  the  shock-on-lip  condition 
and  operating  in  the  self-powered  mode. 


4.1.  The  model 

As  in  Refs.  15,  36,  and  37,  we  consider  hypersonic  gas  flow  along  a  series  of  compression 
ramps  upstream  of  the  inlet  with  a  forward-shifted  cowl  lip.  The  flow  is  two-dimensional  in 
(x,  z )  plane.  Cases  both  without  and  with  MHD  influence  on  the  flow  were  computed.  In  MHD 

cases,  both  magnetic  field  and  ionizing  electron  beam  were  directed  parallel  to  z-axis.  Because 
the  entire  flow  region  is  hypersonic,  steady  state  solution  using  x  as  marching  coordinate  can  be 
found. 

The  set  of  steady-state  Euler  equations  in  Cartesian  coordinates,  together  with  an  ideal  gas 
equation  of  state,  a  simple  model  of  an  ideal  Faraday  MHD  generator,  vibrational  relaxation,  and 
a  plasma  kinetic  model  are  those  of  Refs.  15,  36,  and  37.  Two  significant  changes  were  made  to 
the  code  used  in  Ref.  36.  First,  the  potentiality  of  electric  field,  Vx£  =  o,  was  enforced,  which, 
combined  with  nonuniformity  of  flow  parameters  and  magnetic  field,  results  in  non-trivial 
corrections.  Indeed,  consider  a  flow  at  a  velocity  u  along  x  coordinate.  The  magnetic  field  is 
B=BZ,  and  the  electrical  conductivity  is  <x  In  general,  u  =  u(x,z);  B  =  B (x, z);  a  =  a(x,z) .  The  width 


of  MHD  region  along  y  will  be  denoted  as  H. 

The  Faraday  current  density  is: 
jv  (x,z)  =  a(x,z)(u(x,z)B(x,z)-Ey)  (14) 

where  the  electric  field  £/=const,  independent  of x  and z,  to  satisfy  the  condition  Vxl  =  o. 
The  total  current  collected  on  an  external  resistor  R  is: 


I  =  ||  jv  (x,  z)  dxdz 

x,z 

=  jj  cr(x,z')(u(x,z')B(x,z')  —  Ey'jdxdz  (15) 

X,Z 

=  ||  auBdxdz  -  Ev  ||  a  dxdz 

x,z  x,z 

The  Ohm’s  law  for  the  closed  circuit  can  be  written  as: 


u(x,z)B(x,z)H 


Jy(x’z ) 
a  (x,z) 


H  +  IR  . 


(16) 


Substituting  jy(x,z )  from  (14)  and  /  from  (15)  into  (16),  and  solving  the  resulting  equation  for 
Ey,  we  obtain: 

7?||  auBdxdz 

Er  =  x’z  - . 

H  +  M|  a  dxdz 


(17) 
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The  total  electrical  resistance  of  the  nonunifonn  plasma  region  is: 
H 


||  adxdz 


(18) 


and  the  load  factor  is: 

,  R  R 

k  = 


R  +  r  H 


(19) 


||  adxdz 


From  (19),  the  external  resistance  R  can  be  expressed  in  tenns  of  the  plasma  conductivity  and  the 
load  factor: 


R  =  —H~n-^ - . 

1  -k  1 1  adxdz 


(20) 


Inserting  (20)  into  (17),  we  obtain  the  formula  for  electric  field: 

auBdxdz 


E.  =- 


||  adxdz 


(21) 


This  fonnula  is  a  generalization  for  non-uniform  cases  of  the  familiar  expression  Ey  =  kuB  used 

in  the  case  of  uniform  flow  properties  and  magnetic  fields. 

MHD  computations  start  with  calculation  of  the  electric  field  with  formula  (9).  Then,  local 
current  density  is  be  found  from  (14).  Since  the  conductivity  and  velocity  profiles  in  the  MHD 
region  can  only  be  detennined  from  the  computations,  the  computational  procedure  is,  by 
necessity,  iterative. 

Since  the  field  (21)  is  essentially  an  average  value  of  uB,  multiplied  by  k,  then,  according 
to  Eq.  (14),  in  the  regions  where  the  local  uB  is  less  than  the  average  value  times  k,  current  will 
flow  in  the  reverse  direction. 

External  parameters  of  the  MHD  circuit,  the  current  /  and  the  voltage  V,  are  calculated 
easily: 

/  =  (!-£)  ||  auBdxdz  (22) 


V  =  E  ,H  =  H- 

y 


£||  auBdxdz 

x,z 

||  adxdz 


(23) 


The  second  modification  of  the  computational  procedure  relates  to  modeling  of  power 

15  28  30  50  53 

deposition  and  ionization  by  the  electron  beams.  In  our  earlier  work,  '  ’  '  we  used  the  so- 
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called  ‘forward-backward’  approximation54  for  electron  beam  propagation  along  magnetic  field 
lines  into  gases  and  beam-generated  ionization.  More  recently,  having  analyzed  the  physics  of 
ionization  processes  and  results  of  experiments  and  computations,  we  concluded  that  in  a 
uniform  gas  the  ionization  rate  profile  is  to  be  close  to  a  truncated  Gaussian.  ’  The  simplified 
model  was  then  developed  and  used  for  MHD  analyses.  ’  However,  the  Gaussian  model 
assumed  that  the  electron  beam  is  monoenergetic  at  its  injection  point.  If  the  injected  beam 
electrons  have  a  range  of  energies,  the  resulting  power  deposition  and  ionization  profiles  can  be 
obtained  as  linear  superpositions  of  the  Gaussian  profiles  with  appropriate  weights.  To  turn  this 
statement  around,  any  beam  power  deposition  profile  can  be  represented  as  a  linear  combination 
of  a  finite  or  infinite  number  of  Gaussian  profiles,  corresponding  to  a  certain  initial  energy 
spectrum  of  the  beam  electrons. 

Thus,  we  can,  in  principle,  select  a  beam  power  deposition  profile  that  optimizes  the  MHD 
perfonnance,  which  is  done  in  this  paper,  and  only  then  calculate  the  required  beam  energy 
spectrum  at  the  injection  point  (this  part  is  left  for  a  separate  future  study).  In  fact,  our 
preliminary  analysis  showed  that  the  best  performance  of  on-ramp  MHD  flow/shock  control 
devices  is  reached  with  a  near-uniform  beam  power  deposition  profiles  in  the  vertical  (x) 
direction.  Indeed,  if  the  beam-induced  ionization  profile  peaks  close  to  the  ramp  surface,  then 
little  MHD  effect  exists  far  from  the  surface,  resulting  in  little  effect  on  the  bow  shock  angle  and 
cowl  incidence  point.  If,  on  the  other  hand,  ionization  peaks  far  from  the  surface,  the  effect  on 
the  shock  is  substantial,  but  since  there  is  little  ionization  in  the  region  of  strong  magnetic  field 
(near  the  surface),  the  extracted  MHD  power  is  reduced  and  can  be  insufficient  to  run  the 
ionizing  electron  beams.  Therefore,  in  all  calculations  presented  in  this  paper,  we  used  a  uniform 
profile  of  electron  beam  power  deposition  throughout  the  MHD  region,  characterizing  the  beams 
with  power  deposition  density,  Qb,  and  dimensions  of  the  ionized  region  along  x  and  z. 

For  perfonnance  assessment  of  the  propulsion  system  with  and  without  MHD  control  (the 
detailed  assessment  is  outside  the  scope  of  the  present  work),  the  following  dimensionless 
parameters  were  computed  at  the  inlet  throat:  mass  capture  ratio  km  ;  kp  ;  total  enthalpy  ratio  kH  ; 

static  pressure,  density,  and  temperature  ratios  kp,  kp ,  and  k, ;  cooled  kinetic  energy  efficiency, 
71KE  cooi  ;  and  the  average  Mach  number  M .  The  coefficients  k  for  mass  flow  rate  and  static 
pressures,  temperatures,  etc.,  are  the  ratios  of  the  respective  dimensional  values  at  the  throat  to 
the  freestream  values.  Note  that  mass  capture  km  is  calculated  by  referencing  the  actual  mass 
flow  rate  at  the  throat  to  the  mass  flow  rate  through  the  inlet  capture  area  at  zero  angle-of-attack. 
With  this  definition,  km  can  exceed  1 . 
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In  most  cases,  the  flow  at  the  inlet  throat  is  quite  non-unifonn.  In  this  work,  we  used  the 
so-called  stream-thrust  averaging  commonly  accepted  in  inlet  design.14  This  procedure, 
described  in  Ref.  14,  effectively  takes  into  account  losses  of  total  pressure  (entropy  increase)  that 
would  occur  in  the  isolator  when  the  flow  is  allowed  to  settle  and  to  become  uniform  and  parallel 
to  the  walls. 

4.2.  MHD  control  at  M>MriP<jPn:  methodology  of  computations,  computed  cases, 

and  results 

The  two-dimensional  four-ramp  inlet  geometry,  ’  designed  for  Mach  5  flight  at  2-degree 
angle-of-attack,  had  the  ramp  angles  of  2.5,  8.5,  11,  and  13  degrees.  The  location  of  the  cowl  lip 
was  chosen  so  that  the  first  three  oblique  shocks  would  together  reach  the  point  slightly  upstream 
of  the  lip  (barely  missing  the  lip);  the  fourth  shock  was  allowed  to  hit  the  cowl.  The  freestream 
conditions  at  both  design  and  off-design  (Mach  8)  Mach  numbers  studied  in  this  paper 
correspond  to  the  flight  dynamic  pressure  of  1000  psf  (about  0.5  atm).  Table  1  lists  freestream 
conditions:  altitude  h,  and  static  pressure  and  temperature,  p0  and  T„ ,  in  all  computed  cases. 

At  Mach  numbers  greater  than  5,  the  shocks  would  make  contact  with  the  cowl,  reducing 
total  pressure  at  the  inlet,  and  potentially  causing  very  high  heat  transfer  rates  at  the  shock 
impingement  point,  and  possibly  flow  separation  and  engine  unstart. 

To  maximize  MHD  control  performance,  the  magnetic  field  should  be  as  strong  as 
possible,  and  it  should  protrude  from  the  ramp  into  the  flow  as  far  as  possible.  The  field  was 
assumed  to  be  generated  by  a  superconducting  coil  placed  inside  the  forebody  and  projecting 
magnetic  field  downward.  The  field  strength  at  the  ramp  surface  was  assumed  to  be  3-5  Tesla. 
The  field  outside  the  coil  is  a  function  of  the  ratio  of  the  distance  from  the  coil  end  plane  and  the 
coil  radius.  Thus,  the  protrusion  of  the  B  field  into  the  gas  increases  with  the  coil  diameter. 
Obviously,  to  make  the  magnetic  field  reasonably  unifonn  in  the  spanwise  (y)  direction,  non¬ 
circular  magnets  are  needed.  In  the  present  two-dimensional  modeling,  however,  we  used  the  B 
field  strength  as  a  circular-coil  function  of  z,  with  the  ‘effective’  coil  radius  R  that  should  be 
minimized  in  the  modeling. 

Off  the  centerline,  the  B  field  is  not  only  reduced  in  magnitude,  but  it  also  diverges. 
Therefore,  to  maximize  the  ponderomotive  forces  pushing  the  flow  and  shocks  upstream,  the 
MHD  interaction  region  should  be  concentrated  near  the  centerline  of  the  coil. 

The  ionization  of  air  in  the  MHD  region  is  created  by  electron  beams  injected  from  the 
ramp  along  magnetic  field  lines.  Concentrating  electron  beam  current  in  a  short  region  near  the 
magnet  centerline  would  avoid  the  need  to  have  a  large  area  of  the  vehicle  surface  covered  by 
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fragile  windows  or  differentially  pumped  ports,  which  would  have  been  required  in  cases  of 
long,  distributed  MHD  interaction  regions.  As  the  MHD  interaction  length  L  becomes  short,  the 
electrical  conductivity  <r  has  to  be  increased  in  order  to  keep  the  interaction  parameter  (Stuart 
number)  at  the  same  level  as  that  for  a  long,  distributed  MHD  region.  Since  the  ionization 
fraction  and  the  conductivity  are  proportional  to  the  square  root  of  electron  beam  current  density, 
shortening  the  MHD  region  from  2-3  meters  down  to  10-20  cm  requires  increasing  the  beam 
current  density  by  2  orders  of  magnitude,  to  50-100  mA/cm  .  Aerodynamic  windows  or 
differentially  pumped  plasma  portholes  may  be  good  candidates  for  electron  beam  transmission. 
Such  windows  or  ports  can  presumably  handle  very  high  electron  beam  currents,  even  much 
higher  than  100  mA/cm",  in  contrast  to  conventional  thin  foils. 

As  discussed  in  Section  4.1,  the  beam  power  deposition  profde  Qh  was  assumed  to  be 
unifonn  in  the  entire  MHD  region:  xMHD  <  x  <  xMHD  +  AxMHD ;  \z-zb(x)\  <  Lb . 

The  goal  of  optimization  studies  was  to  find  the  values  of  parameters  corresponding  to 
restoration  of  the  design  shock  wave  configuration  at  minimal  magnet  radius,  with  MHD 
extracted  power  satisfying  the  constraint  Pb  ~\PMHD  ■  The  latter  criterion  expresses  the  condition 

that  the  MHD  device  be  self-powered,  i.e.  that  the  MHD  generated  power  be  just  enough  to 
generate  the  ionizing  beams,  given  the  approximately  50%  efficiency  of  electron  beam 
generation. 

At  a  given  MHD  generator  position  along  the  forebody  and  B-field  strength  at  the  ramp 
surface,  optimization  parameters  are:  effective  radius  of  the  magnet,  R;  e-beam  penetration 
depth,  Lh ;  and  MHD  channel  length,  Axv,„„ . 

For  the  best  perfonnance,  the  MHD  interaction  region  should  be  placed  as  far  upstream  as 
possible.  ’  Indeed,  the  required  linear  shift  of  the  shock  incidence  point  at  the  cowl  level  is 
caused  by  a  smaller  change  in  shock  angle,  and,  therefore,  by  a  lower  MHD  interaction 
parameter,  when  the  MHD  region  is  moved  upstream.  Additionally,  creating  the  necessary  level 
of  ionization  at  lower  gas  density  (closer  to  freestream  conditions)  requires  lower  current  density 
and  power  of  the  ionizing  electron  beams.  The  two  requirements  for  the  placement  of  the  MHD 
interaction  region  (that  the  interaction  region  should  be  as  far  upstream  as  possible,  and  that  it 
should  be  close  to  the  magnetic  coil  centerline)  define  the  minimum  distance  of  about  one 
magnet  radius  of  the  MHD  interaction  region  from  the  nose. 

Thus,  Cases  A  (Ai-A4),  B,  and  C  correspond  to  MHD  region  placed  on  the  first  ramp,  at 
the  minimal  distance  from  the  nose.  The  principal  parameters  in  these  and  other  cases  are  listed 
in  Table  2.  The  inlet  performance  predictions  are  listed  in  Table  4.  Table  3  contains  the  principal 
MHD  parameters:  the  rate  of  work  by  j  x  B  forces,  P.  s,  the  extracted  power,  PMHD,  the  gas 
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heating  rate,  P, ,  the  power  deposition  into  vibrational  mode  of  molecules,  Pv ,  the  power  of 
ionizing  electron  beams,  Ph ,  the  MHD  interaction  parameter,  S  ,  the  enthalpy  extraction  ratio,  77 , 
the  total  Faraday  current,  I  ,  and  the  Faraday  voltage,  Uel .  Due  to  the  two-dimensionality  of  the 
problem,  all  power  quantities  P  and  Uel  are  expressed  per  unit  length  ( 1  meter)  in  the  spanwise 
(y)  direction.  Note  that  the  total  rate  of  work  by  j  x  b  forces,  P]xS,  is  spent  on  1)  irreversible 
Joule  heating  that  consists  of  gas  heating  in  the  narrow  sense,  P, ,  and  the  power  deposition  into 
vibrational  mode  of  molecules,  Pv ,  and  2)  the  generated  electric  power,  PMHD : 
P]xS=PMHD+(Pj+Pv).  The  ratio  between  the  extracted  power  and  the  rate  of  dissipation  is 
determined  by  the  load  factor.  At  k=0.5,  PMHD  =  P,  +  Pv  =  0.5  P]xS .  Note  also  that  PMHD  =  ly  x UeI . 

To  demonstrate  advantages  of  upstream  placement  of  the  MHD  region,  Cases  D  and  E 
were  computed  with  MHD  region  shifted  downstream.  As  seen  in  Table  2,  both  R  and  Lb  have  to 
be  increased  in  these  cases  compared  with  those  in  Cases  Ai  and  C.  Additionally,  due  to  the 
more  intense  flow  heating  and  the  need  to  have  a  larger  e-beam  ionized  region  (longer  Lb),  flow 
spillage  and  reduction  in  air  mass  capture  occur  in  Cases  D  and  E  (see  Table  4).  Also,  although 
Cases  D  and  E  yield  higher  compression  ratios  kp  than  Cases  Ai  and  C,  the  better  compression 
comes  at  the  expense  of  reduced  total  pressure  and  riKi:  cool .  Note  that  our  numerous  calculations 

with  variation  of  all  parameters  showed  that  for  the  given  design  and  flight  conditions,  it  is 
impossible  to  restore  the  shock-on-lip  condition  at  Mach  8  with  MHD  region  located  at  x>13  m, 
even  with  a  ‘softened’  self-power  condition  PMHD  >  Pb . 

Cases  Ai,  B,  and  C  clearly  demonstrate  the  trade-off  between  the  strength  of  magnetic  field 
and  the  size  of  the  magnet.  Comparison  of  these  three  cases  shows  that  the  magnet  size  can  be 
reduced  and  the  shock-on-lip  condition  restored,  if  the  strength  of  magnetic  field  is  increased 
(Table  2).  However,  as  seen  in  Table  4,  Case  Ai  is  the  best  among  these  three  cases  in  terms  of 
compression  ratio  kp,  kinetic  energy  efficiency  rjKEcool ,  and  enthalpy  ratio  kH.  Thus,  larger 

magnets  with  relatively  modest  B  field  strength  do  work  better  than  smaller  magnets  with 
stronger  B. 

All  cases  except  A2  had  load  factors  k=0.5.  Case  A2  was  run  at  a  higher  load  factor, 
k=0.75.  As  seen  in  Table  2,  increasing  load  factor  requires  larger  magnets  and  longer  MHD 
regions.  Calculations  at  k  less  than  0.5,  not  listed  in  the  tables,  showed  that  the  performance 
becomes  worse  than  that  at  k=0.5,  because  with  less  power  extraction  it  is  difficult  to  restore  the 
design  shock  configuration,  while  demanding  that  Pb  » \PMHD  .  Thus,  for  optimal  performance,  the 

load  factor  should  not  deviate  far  from  0.5. 
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Comparison  of  Cases  Ai,  A3,  and  A4  demonstrates  the  role  of  the  magnitude  of  Qh 
(electron  beam  power  deposition  density).  As  seen  in  Table  2,  at  the  given  flight  conditions,  the 
value  Qb  =  10  MW/m  (Case  Ai)  is  close  enough  to  the  optimum  for  the  MHD  region  to  be  short, 
Axmhd  <  0.3  m,  while  satisfying  the  criterion  PMHD  » \Ph  .At  Qb=2  MW/nr  (case  A3),  shock-on-lip 
condition  can  be  achieved  at  substantially  increased  R,  Lb,  and  AxMHD .  At  Qb  =50  MW/m  (case 
A4),  the  values  of  R,  Lb  and  AxMHD  can  be  decreased  as  compared  with  those  in  the  Case  Ai,  but 
only  if  the  criterion  PMHD  ~\Pb  is  relaxed  to  PMHD  ~  Ph . 

Figures  7  and  8  show  the  location  of  the  MFID  region,  contour  lines  of  gas  and  vibrational 
temperature,  and  flow  streamlines  in  Cases  Ai  and  D.  As  seen  in  the  figures,  due  to  vibrational 
excitation  by  electron  impact  in  the  MHD  region  and  relatively  slow  vibrational  relaxation  in  the 
cold,  low-density,  weakly  dissociated  flow,  there  is  a  considerable  vibrational  nonequilibrium  at 
the  inlet  throat.  One  challenge  posed  by  this  nonequilibrium  is  that  the  stream-thrust  averaging 
procedure14  was  developed  for  thermally  equilibrium  flows.  Fortunately,  the  kinetic  energy 
efficiency  nKE,cool  was  found  to  be  not  very  sensitive  to  the  nonequilibrium  conditions. 

Nevertheless,  the  procedure  of  stream-  thrust  averaging  should  be  generalized  for  nonequilibrium 
flow  in  the  future. 

We  assumed  the  MHD  region  to  be  an  ideal  Faraday  generator.  The  latter  assumption  is 
non-trivial.  The  Hall  parameters  in  the  computed  cases  can  reach  as  high  as  10.  In  conventional 
Faraday  channels  with  continuous  electrodes,  MHD  effects  are  sharply  reduced  due  to  the  Hall 
effect,  while  in  Faraday  channels  with  segmented  electrodes  there  is  a  danger  of  arcing  between 
the  electrode  segments  at  high  Hall  parameters.49  However,  the  MHD  devices  analyzed  in  this 
work  are  very  different  from  conventional  devices.  Not  only  is  the  conductivity  non-thennal, 
generated  and  controlled  by  externally  injected  electron  beams,  but  also  the  length  of  the  channel 
(15-30  cm)  is  much  shorter  than  its  width  (5  meters).  Considering  short  MHD  channel  as  a 
“slice”  of  a  long  channel,  with  only  1  electrode  segment,  one  might  argue  that  the  performance 
of  such  a  short  channel  should  be  close  to  that  of  an  ideal  Faraday  device,  even  at  high  Hall 
parameters.  This  argument  requires  an  investigation  and  validation. 

Overall,  MHD  control  of  hypersonic  flows  appears  promising.  Note  that  the  first 
experimental  demonstration  of  cold-air  supersonic  MHD  effect  with  non-self-sustained 
ionization  (by  repetitive  high-voltage  nanosecond  pulses)  has  been  recently  performed  at 
Princeton  University.55 
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Table  1.  Freestream  conditions  in  computed  cases.  Freestream  dynamic  pressure  is 

q=1000  psf. 


Mach  number 

h,  km 

P0,  Pa 

To,  K 

5 

24.538 

2735 

221.09 

8 

30.76 

1068.75 

227.26 

Table  2.  MHD  generator 
cases 


Cas 

Qb, 

B0,  T 

R,  m 

L„,  m 

XMHD 

Ar  MHD  ,  lai 

k 

e 

MW/m3 

A, 

10 

3 

1.65 

1.65 

2 

0.3 

0.5 

a2 

10 

3 

1.75 

1.65 

2 

0.375 

0.75 

a3 

2 

3 

2 

2.3 

2 

0.65 

0.5 

a4 

50 

3 

1.3 

1 

2 

0.15 

0.5 

B 

10 

4 

1.15 

1.5 

1.25 

0.25 

0.5 

C 

10 

5 

0.8 

1.25 

1 

0.25 

0.5 

D 

10 

3 

2.6 

2.5 

9 

0.25 

0.5 

E 

10 

5 

1.5 

2.25 

9 

0.25 

0.5 

Table3. Computed  MHD 
parameters. 


Case 

PhS 

MW/m 

Pmhd  M 

W/m 

Pj, 

MW/m 

Py, 

MW/m 

P„, 

MW/m 

5 

n 

V 

kA 

Uel, 

kV/m 

A! 

25.4 

10 

9.21 

6.218 

4.926 

0.12 

0.02 

6.66 

1.5 

a2 

22.8 

11.91 

7.53 

3.41 

6.14 

0.19 

0.024 

5.146 

2.315 

a3 

32.6 

11.89 

1.19 

8.78 

2.97 

0.1 

0.024 

8.67 

1.37 

a4 

7.36 

7.94 

6.62 

3.68 

7.35 

0.15 

0.016 

4.6 

1.72 

B 

23.33 

7.82 

8.33 

7.17 

3.72 

0.19 

0.016 

5.23 

1.49 

C 

20.83 

6.07 

7.74 

7.014 

3.098 

0.12 

0.012 

4.27 

1.42 

D 

36.71 

14.88 

15.03 

6.8 

7.47 

0.1 

0.03 

9.52 

1.56 

E 

46.1 

14.03 

18.29 

13.76 

6.72 

0.12 

0.28 

9.0 

1.55 
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Table  4.  Computed  inlet  parameters  with  and  without  MHD  control  for  Mach  5  design. 
Angle  of  attack  is  2°;  flight  dynamic  pressure  is  q=1000  psf  for  all  MHD  cases. 


Freestream 

Mach  No. 

and  MHD 

conditions 

K 

kf 

kP 

kT 

M 

KE, cool 

kH 

: 

de: 

1 

M=5; 
sign,  no 
V1HD 

1.14 

16.07 

6.94 

2.31 

2.822 

0.991 

1 

8 

No 

MHD 

1.13 

22.01 

6.547 

3.36 

3.94 

0.969 

1 

Ai 

1.14 

22.46 

6.8 

3.3 

3.87 

0.924 

0.955 

a2 

1.14 

21.92 

6.78 

3.23 

3.92 

0.926 

0.956 

a3 

1.13 

4 

23.09 

6.83 

3.37 

3.78 

0.911 

0.945 

A4 

1.14 

23.18 

6.75 

3.43 

3.82 

0.943 

0.976 

B 

1.14 

20.39 

6.773 

3.01 

4.07 

0.919 

0.944 

C 

1.14 

19.66 

6.74 

2.91 

4.16 

0.92 

0.943 

D 

1.09 

30.41 

6.85 

4.43 

3.16 

0.91 

0.966 

E 

1.11 

35.62 

7.26 

4.9 

2.89 

0.886 

0.95 

4.3. The  effect  of  current  reversal 

As  discussed  in  Section  4.1,  in  those  areas  within  the  MHD  region  where  the  local 
value  of  uB  is  less  than  the  average  value  multiplied  by  k,  the  electric  current  should  flow 
in  the  direction  opposite  to  the  Faraday  electromotive  force.  This  effect  was  indeed 
observed  in  all  computed  cases.  As  an  example,  Fig.  9  shows  the  current  profile  in  Case 
E.  As  seen  in  the  Figure,  the  current  density  reduces  to  zero  far  from  the  ramp,  and  then 
becomes  negative.  In  the  region  of  ‘negative’  current,  the  jxB  force  reverses  direction 
and  becomes  an  accelerating  force.  Although  the  effect  was  seen  in  all  computed  cases, 
the  results  were  not  changed  much  by  it.  However,  it  is  conceivable  that  in  some  other 
MHD  problems  the  current  reversal  can  play  a  significant  role. 

Consider,  for  example,  the  near-wall  region  of  an  MHD  generator.  If  the  B  field  is 
nonnal  to  the  surface,  and  the  conductivity  is  more  or  less  uniform  across  the  boundary 
layer,  then  close  enough  to  the  wall,  where  the  velocity  is  less  than  the  freestream 
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velocity  multiplied  by  the  load  factor,  the  jxB  will  accelerate  the  boundary  layer.  This 
effect  can,  in  principle,  change  the  boundary  layer  structure  and  thickness,  and  can  affect 
vorticity  generation,  flow  separation,  and  laminar-turbulent  transition.  These  phenomena 
should  be  investigated  in  future. 


] 


Fig.  7.  Location  of  the  MHD  interaction  region,  contours  of  gas  and  vibrational 
temperatures,  and  streamlines  at  Mach  8  with  MHD  control  for  the  case  of  magnet  with 
B0=3  T;  R=1.65  m  (Case  Ai) 
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Fig.  8.  Location  of  the  MHD  interaction  region,  contours  of  gas  and  vibrational 
temperatures,  and  streamlines  at  Mach  8  with  MHD  control  for  the  case  of  magnet  with 
B0=3  T;  R=2.6  m  (Case  D) 
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Fig.  9.  Current  density  profile  and  directions  in  the  MHD  region  for  Case  E. 
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5.  Air  capture  increase  with  Virtual  Cowl 

The  Virtual  Cowl  concept  15, 35  is  schematically  shown  in  Fig.  10.  We  considered 
2D  hypersonic  gas  flow  along  a  series  of  compression  ramps  upstream  of  the  inlet;  cases 
both  without  and  with  heat  addition  were  computed  with  the  set  of  Euler  equations 
together  with  perfect  gas  equation  of  state.  Because  the  entire  flow  region  is  supersonic,  a 
steady  state  solution  using  x  as  marching  coordinate  can  be  found.  The  heat  addition  rate 
profile  was  set  as  Gaussian,  that  is,  it  was  proportional  to 

exp(-(jc-x0)2/^,x -{z-zaf  ^  ,  where  x0  and  z0  are  coordinates  of  the  center  of  the 
heating  region,  and  reff  x  and  reff  __  are  the  effective  radii.  The  parameters  x0 ,  z0 ,  reff  x ,  and 
r  „  were  varied  in  computations  in  order  to  find  optimal  values  of  these  parameters.  In 
“circular-cylinder  Gaussian”  cases,  reff  x  =  reff  z  =  reff .  In  “elliptical-cylinder  Gaussian” 
cases,  reff  x  *  reff  z ,  the  heating  profile  may  also  be  tilted  at  an  angle  aq  with  respect  to  x 
axis. 

The  inlet  performance  was  characterized  by  the  set  of  stream-thrust  averaged 
parameters,  as  in  the  MHD  control  studies  in  Section  4. 


Fig.  10.  Schematic  diagram  of  the  virtual  cowl  concept:  off-body  heat  addition  increases 
mass  capture. 

Freestream  conditions  in  computed  cases  corresponded  to  flight  at  Mach  6,  8,  or  10, 
with  dynamic  pressure  of  either  1000  psf  or  2000  psf.  The  baseline  case  was  that  of  Mach 
10  flight  with  zero-degree  angle-of-attack.  The  3-ramp  inlet  (2.5°,  8.0°,  and  11.0°  ramp 
angles)  was  chosen  so  that  the  nose  shock  would  barely  (by  about  1  inch)  miss  the  cowl 
lip,  and  the  flow  would  be  almost  parallel  to  the  isolator  walls  downstream  of  the  inlet 
throat. 
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To  reduce  the  spillage,  the  cowl  might  be  extended  upstream,  although  this  could 
be  difficult  to  accomplish  practically.  However,  the  case  with  extended  solid  cowl  can 
serve  as  a  guide  to  optimizing  virtual  cowl  cases.  Qualitatively,  one  would  expect  that  the 
virtual  cowl  (i.e.,  off-body  energy  addition)  should,  for  the  best  performance,  create  a 
flow  pattern  imitating  that  of  the  solid  extended  cowl.  Since  shocks  do  not  reflect  off  a 
heated  region,  the  analogy  between  the  solid  and  virtual  cowls  is  not  exact.  However,  the 
heated  region  can  generate  a  shock  wave  that  would  mimic  the  shock  reflected  from  the 
solid  surface  upon  incidence  of  the  nose  shock  on  it.  Thus,  one  would  expect  that  the 
energy  addition  should  be  more  or  less  concentrated,  and  that  it  should  be  positioned  near 
the  intersection  of  the  nose  shock  with  the  upstream  continuation  of  the  cowl. 

Note  that  heating  and  the  resulting  expansion  of  the  gas  would  adversely  affect  both 
mass  flow  rate  and  stagnation  pressure  at  the  throat  if  the  heated  air  is  allowed  to  enter 
the  inlet.  Therefore,  the  energy  addition  region  should  be  shifted  somewhat  down,  so  that 
heated  and  expanded  air  would  miss  the  inlet,  while  cold  air  would  be  compressed  and 
deflected  into  the  inlet. 

The  qualitative  expectations  of  the  preceding  two  paragraphs  were  confirmed  in 
calculations.  With  “circular-cylinder  Gaussian”  heat  addition,  reff  x  =  reff  z  =  reff  =  0.2  m  was 

found  to  be  close  to  the  optimum.  Calculations  showed  that  energy  addition  does  increase 
the  mass  capture,  and  the  effect  becomes  stronger  as  more  energy  is  added.  With  varying 
coordinate  of  the  heated  region,  the  mass  capture  and  the  enthalpy  flux  both  exhibit  broad 
maxima  at  the  value  of  x  approximately  corresponding  to  the  intersection  of  the  nose 
shock  with  the  upstream  continuation  of  the  cowl.  Total  pressure  and  adiabatic  kinetic 
energy  coefficients  also  have  their  peaks,  but  at  a  value  of  x  somewhat  shifted 
downstream  with  respect  to  that  where  mass  capture  is  maximum.  As  predicted  in  the 
previous  paragraph,  shifting  the  center  of  energy  addition  region  by  2 reff  =  0.4  m  down 

from  the  cowl  line  somewhat  improves  inlet  performance.  Static  temperature  contours 
and  flow  streamlines  in  Mach  6,  1000  psf  case  with  energy  addition  rate  of  W= TO  MW/m 
centered  near  the  optimum  location  are  shown  in  Fig.  11. 
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Fig.  11.  Static  temperature  contour  lines  (upper  plot)  and  flow  streamlines  (lower  plot)  at 
Mach  6,  1000  psf,  with  optimum  location  of  “circular  Gaussian”  Virtual  Cowl  with 
reff  =  0.2  m  ,  W=  10  MW/m. 


We  found  that  if,  instead  of  a  single  energy  addition  region,  two  energy  addition 
regions  are  used  with  the  total  power  equal  to  that  of  the  single  region,  then,  with  proper 
positioning  of  the  second  source  downstream  of  the  first  one,  the  inlet  performance  would 

"3  c 

not  decrease  compared  with  that  given  by  the  single  source. 

We  have  also  run  a  series  of  calculations  with  “elliptical-cylinder  Gaussian” 
heating  profiles,  reff  x  *  re// . .  Stretching  the  energy  addition  profile  along  the  flow,  while 

maintaining  constant  total  power,  was  found  to  result  in  a  small  improvement  in  inlet 
performance.  Specifically,  reff  x  =  15 reffz  was  found  to  give  the  best  perfonnance.  We 

further  found  that  slightly  tilting  the  stretched  heating  profile  with  respect  to  x  axis  can 
improve  the  performance.  The  optimum  angle  turned  out  to  be  aq  =  2.5° ,  which,  perhaps 

not  coincidentally,  equals  the  angle  of  the  first  ramp.  Fig.  12  illustrates  the  optimum  case 
with  stretched  and  tilted  heating  profile  at  Mach  6.  Overall,  although  stretched  and  tilted 
heating  profiles  do  perform  better  than  the  “circular”  ones,  the  improvement  is  only 
incremental.35 
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Fig.  12.  Upper  plot:  power  density  of  the  optimally  located,  stretched,  and  tilted  10 
MW/m  source.  Middle  and  lower  plot:  static  temperature  contours  and  flow  streamlines 
at  Mach  6,  1000  psf,  with  that  heating  source. 

Finally,  the  constraint  that  the  energy  addition  has  to  be  placed  at  or  only  slightly 

-j  ^ 

below  the  cowl  level  was  relaxed.  For  the  Mach  6,  1000  psf  flight,  the  “circular- 
cylinder  Gaussian”  heating  profde  with  reff  =  0.2  m  and  the  total  power  of  5  MW/m  was 

positioned  at  different  vertical  locations  z  below  the  cowl  level  zc\.  The  streamwise 
coordinate  of  the  center  of  energy  addition  region  was  held  constant  at 
xQ  =  0.5,  5.5,  and  10.5  m .  In  each  case,  the  best  inlet  perfonnance  turned  out  to  occur  when 
the  shock  generated  by  the  energy  addition  passed  through  the  point  of  intersection  of  the 
nose  shock  with  the  upstream  continuation  of  the  cowl  line,  as  expected.  The  best  ijKi: 
achieved  was  equal,  and  the  best  km  only  incrementally  superior  to  those  in  the  case  with 
energy  addition  slightly  below  the  cowl  level.  Given  that  the  performance  increase  is 
small,  and  that  it  can  be  practically  difficult  to  implement  and  control  energy  addition  far 
from  the  vehicle  surface,  it  appear  advantageous  to  work  with  energy  addition  at  or 
slightly  below  the  cowl  level.  It  is  instructive,  however,  to  note  that,  although  matching 
the  krn  given  by  the  extended  solid  cowl  at  Mach  6  does  not  seem  possible  with  energy 
addition  concentrated  at  or  just  below  the  cowl  line,  increasing  the  power  of  the  heating 
source  while  moving  the  source  away  from  the  body  can  yield  the  km  equal  to  that  of  the 
solid  cowl.  Indeed,  having  run  cases  with  power  depositions  of  20  and  35  MW/m,  we 
found  that  the  35  MW/m  case  does  match  the  km  of  extended  solid  cowl,  although  with 
somewhat  lower  rjKE.  Note  that  at  power  levels  as  high  as  35  MW/m,  the  optimum 
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location  of  the  energy  addition  should  be  farther  from  the  body;  this  was  not  fully 
computed  in  the  present  work. 

The  energy  required  for  the  Virtual  Cowl  may  be  supplied  by  an  in-stream  MHD 
power  extraction  system,  in  which  case  the  additional  mass  flow  may  come  at  the 
expense  of  enthalpy  taken  from  the  propulsive  flowpath.  In  this  work,  sample  scramjet 
engine  cycles  were  used  to  assess  the  impact  of  additional  captured  flow  and  enthalpy 
extraction. 

The  following  six  cases  were  analyzed.  Cases  A  and  B  represent  two  Mach  8 
conditions  with  two  levels  of  optimized  “circular  Gaussian”  energy  deposition,  5  MW/m 
and  10  MW/m,  respectively.  Cases  C  to  F  represent  Mach  6  conditions  at  different 
energy  deposition  levels  and  locations:  Case  C  corresponds  to  the  optimized  “circular 
Gaussian”  energy  deposition  of  5  MW/m;  Case  D  corresponds  to  the  optimized  “circular 
Gaussian”  energy  deposition  of  10  MW/m;  Case  E  corresponds  to  the  optimized  “tilted 
elliptical  Gaussian”  ( refLx  =  3  m,  re// .  =  0.2  m,  aq  =  2.5°)  energy  deposition  of  10  MW/m;  and 

Case  F  corresponds  to  the  “circular  Gaussian”  energy  deposition  of  35  MW/m  located  far 
from  the  cowl  lip  (xq=18.5  m,  zcl-zq=1  .3  m).  The  results  of  the  inlet  calculations  are 
shown  in  third  through  fifth  columns  of  Table  5.  The  results  show  that  for  energy 
extractions  on  the  order  of  1.2%  to  3.4%,  the  mass  flow  through  the  engine  can  be 
increased  by  7. 1%  to  15.5%  with  little  change  in  the  inlet  efficiency. 

Sample  engine  calculations  were  conducted  at  Mach  6  and  8  for  a  baseline  Mach  10 
design  operating  on  hydrogen  fuel.  The  baseline  engine  calculations  assumed  t|ke  =  0.97, 
combustion  efficiency  =  0.9,  and  nozzle  efficiency  =  0.97.  The  assumed  inlet  kinetic 
energy  efficiency  (t|ke  =  0.97)  was  chosen  to  represent  typical  values  for  an  inlet 
including  friction  effects.  The  inlet  efficiency  results  from  the  Virtual  Cowl  analysis 
were  applied  as  a  delta  to  the  assumed  0.97  value. 

The  effects  of  the  virtual  cowl  on  the  engine  performance  were  assessed  by 
increasing  the  mass  flow,  modifying  the  inlet  kinetic  energy  efficiency,  and  removing  the 
required  enthalpy  from  the  scramjet  combustor.  Note  that  this  analysis  assumes  that  the 
energy  deposited  to  create  the  virtual  cowl  is  obtained  from  the  scramjet  combustor  and 
the  efficiency  of  the  energy  extraction  is  100%. 
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The  results  of  the  engine  cycle  calculations  are  shown  in  the  last  two  columns  of 
Table  5,  where  it  can  be  seen  that  the  engine  thrust  coefficient  CY  increases  nearly 
proportionally  to  the  mass  flow  increase.  The  effect  of  the  virtual  cowl  on  the  engine 
specific  impulse  is  seen  to  be  relatively  small  (i.e.  less  than  1%).  Thus,  the  virtual  cowl 
appears  to  represent  a  system  for  augmenting  the  thrust  of  a  scramjet  engine  with  little 
impact  on  overall  cycle  efficiency.  (Apparently,  the  additional  compression  achieved 
with  the  additional  captured  mass  flow  offsets  the  enthalpy  extraction  from  the 
combustor.) 

The  results  from  the  engine  cycle  calculations  show  that  the  application  of  a  virtual 
cowl  merits  additional  investigation.  While  the  required  energy  deposition  levels  are 
high  (5-35  MW  per  meter  of  engine  width),  and  the  weight  impacts  of  power  extraction 
and  energy  deposition  systems  have  not  been  estimated,  the  ability  to  increase  the  engine 
thrust  of  an  engine  for  an  accelerating  system  can  be  significant. 

Two  additional  sets  of  engine  cycle  calculations  were  conducted  to  further  assess 
the  impact  of  the  virtual  cowl  concept.  In  the  first  set  of  calculations,  the  effects  of  the 
virtual  inlet  were  included  in  the  calculation,  but  no  energy  was  extracted  from  the 
propulsion  flowpath.  This  set  of  calculations  applies  in  the  case  when  the  energy  needed 
for  production  of  the  virtual  cowl  is  either  carried  inside  the  vehicle  or  produced  with  an 
on-board  system  that  does  not  impact  the  engine.  In  the  second  set  of  calculations,  the 
virtual  cowl  was  not  used,  but  the  energy  that  would  be  needed  was  assumed  to  be 
applied  to  the  combustor.  This  set  of  calculations  was  conducted  as  a  reference  and 
would  correspond  to  the  case  when  the  additional  energy  is  available  (either  carried  on¬ 
board  or  produced  separately  from  the  propulsive  flowpath)  and  supplied  directly  to  the 
combustor. 

The  results  from  the  two  sets  of  engine  calculations  are  compared  to  the  virtual 
cowl  results  in  Table  6.  The  results  show  that  the  engine  performance  is  better  for  the 
virtual  cowl  when  the  energy  is  not  extracted  from  the  flowpath.  This  is  the  expected 
result  since  energy  extraction  is  equivalent  to  a  reduction  in  fuel  supplied  to  the 
combustor.  The  last  set  of  calculations  shows  the  impact  of  energy  deposition  within  the 
combustor.  In  this  case,  the  calculations  assume  that  the  energy  is  available  on-board  and 
directly  deposited  in  the  combustor.  The  results  show  that  both  the  thrust  coefficient  and 
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specific  impulse  increase  by  small  percentages.  These  results  show  that  the  principal 
advantage  of  the  virtual  cowl  is  the  increase  in  engine  thrust  (primarily  caused  by  the 
increase  in  captured  airflow)  with  little  change  in  specific  impulse  regardless  of  where 
the  energy  for  production  of  the  virtual  cowl  is  obtained. 


Table  5.  Virtual  Cowl  impact  on  engine 

perfonnance 


Case 

Freestream  Mach 

Number 

Enthalpy 

Extracted 

Mass  Flow 

Increase 

A  CT 

A Isp 

A 

8 

1.2% 

7.1% 

+0.002 

+6.8% 

-0.27% 

B 

8 

2.3% 

11.4% 

+0.002 

+  10.5% 

-0.82% 

C 

6 

1.8% 

7.8% 

+0.004 

+8.0% 

+0.20% 

D 

6 

3.4% 

13.9% 

0.0 

+  14.0% 

+0.05% 

E 

6 

3.4% 

15.5% 

+0.003 

+  15.8% 

+0.27% 

F 

6 

10.2% 

34.5% 

-0.003 

+33.6% 

-0.52% 

Table  6.  Comparison  of  virtual  cowl  with  and  without  energy  extraction,  and  effect  of 
energy  deposition  in  combustor. 


Virtual  Cowl  with 

Energy  Extraction 

Virtual  Cowl  without 

Energy  Extraction 

No  Virtual  Cowl  and 

Energy  Addition 

Case 

Freestream 

Mach 

Number 

A  CT 

A Isp 

tj 

<i 

A Isp 

A  CT 

A Isp 

A 

8 

+6.8% 

-0.21% 

+7.6% 

+0.46% 

+0.78% 

+0.78% 

B 

8 

+10.5% 

-0.82% 

+  12.0% 

+0.57% 

+  1.48% 

+  1.48% 

C 

6 

+8.0% 

+0.20% 

+8.7% 

+0.78% 

+0.63% 

+0.63% 

D 

6 

+14.0% 

+0.05% 

+  15.1% 

+  1.10% 

+  1.18% 

+  1.18% 

E 

6 

+  15.8% 

+0.27% 

+  17.0% 

+  1.30% 

+  1.18% 

+  1.18% 

F 

6 

+33.6% 

-0.52% 

+37.2% 

+2.17% 

+3.50% 

+3.50% 
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